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1. Introduction 

The Thermodynamic Bethe Ansatz (TBA) is an efficient tool to determine the finite- 
size spectrum of two-dimensional relativistic integrable models [1]. The relevance of 
the TBA for the AdSs x spectral problem was advocated in [2] where Liischer's 
approach [3] was used to relate exponential corrections to string energy to wrapping 
effects in dual field theory. Its application to a nonrelativistic model requires under- 
standing thermodynamic properties of a so-called mirror theory [4] which is obtained 
from the original model by means of a double- Wick rotation. The mirror model of 
the light-cone AdSs x superstring was studied in detail in [4] where, in particu- 
lar, the asymptotic spectrum was identified, and the mirror form of the Bethe- Yang 
equations of [5] was determined. 

The major step towards realizing the TBA approach is to formulate a so-called 
string hypothesis [6] which classifies the states contributing in the thermodynamic 
limit. TBA equations and the associated Y-system are then readily derived from 
it, see e.g. [7]. This step was made last year in [8] where the mirror Bethe- Yang 
equations [4] were used to formulate a string hypothesis for the AdSs x mirror 
model. This opened a way to derive the corresponding canonical [9, 10, 11] and 
simplified TBA equations [12]. 

The TBA equations combined with a certain analytic continuation procedure 
proposed for relativistic models in [13, 14] can be also used to find energies of ex- 
cited states.^ The procedure was called the contour deformation trick in [28] because 

^For other approaches and appUcations, see e.g. [15]-[27]. 
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it basically reduces to deforming the integration contours in ground-state TBA equa- 
tions while keeping their form untouched. As a result, excited states TBA equations 
differ only by integration contours. For practical applications however one should 
take the integration contours back to their ground-state positions, and this results 
in the appearance of state-dependent driving terms in the TBA equations. 

The contour deformation trick was used in [28] to analyze two-particle states 
from the 51(2) sector. It was shown there that two-particle states are divided into 
infinitely-many classes, each class having its own set of driving terms in the TBA 
equations. For the Konishi-like states the TBA equations of [28] are believed to be 
equivalent to those of [11, 29],^ and our numerics strongly supports this. 

The TBA equations obtained in [28] were used to derive the five-loop anomalous 
dimension of the Konishi operator, and to show numerically [33] that the correspond- 
ing result perfectly agrees with the one recently obtained via the generalized Liischer 
formulae [34]. The analysis in [28, 33] was then extended to prove the agreement 
analytically first for the Konishi operator [35], and then for an arbitrary twist-two 
operator [36] reproducing the results in [37].^ For other applications of the TBA 
approach to the AdS/CFT spectral problem, see [52]- [57]. 

In a parallel development the TBA equations proposed in [11] were used in [29] 
to compute anomalous dimension of the Konishi operator up to a relatively large 
value of the 't Hooft coupling constant A ~ 664. Analyzing the results obtained, the 
authors of [29] found the following fitting function for the Konishi state energy or, 
equivalently, for the conformal dimension of the dual Konishi operator 



This formula allows one to make four remarkable predictions. First, it predicts that 
the coefficient c_i of the leading term in the large A expansion is equal to 2 which 
agrees with the spectrum of string theory in flat space [58] and asymptotic Bethe 
ansatz considerations [59]. Second, it shows the vanishing of the constant term cq. 
Third, it makes a new prediction which disagrees with the semiclassical consideration 
in [60] that the flrst nonvanishing subleading coefficient is also equal to 2. And fourth, 

^The set of integral equations proposed in [11, 29] was named an integral form of the Y-system. 
We find this terminology somewhat misleading because these integral equations were not derived 
from the Y-system conjectured in [30] but were proposed by following the pure TBA approach, that 
is by using the canonical TBA equations [9, 10, 11], the contour deformation trick [13], and the 
large J asymptotic solution of the Y-system [30]. A derivation of the TBA equations from the Y- 
system requires understanding the complicated analytic structure of Y-functions, see [9, 31, 28, 32] 
for some results in this direction. 

^The five-loop computations in [34, 37] were based on [38] where the four-loop anomalous di- 
mension of the Konishi operator was computed and shown to agree with the direct field theory 
result [39, 40] . This constituted a strong test of the AdS/CFT correspondence [41]. For other 
applications of Liischer's approach, see e.g. [42]-[51]. 
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it predicts that up to an overall factor of the Konishi state energy is a series in 
1/\/A. The last prediction in fact agrees with the argumentation in [60] but it 
disagrees with the considerations in [61] where the free fermion model describing the 
su(l|l) sector in the semi-classical approximation was analyzed, and it was shown 
that the strong coupling expansion for short states is in powers of 1/ v^- It is worth 
noting that this simple model indeed predicts that the constant term in the strong 
coupling expansion vanishes. The formulae derived in the framework of the free 
fermion model would definitely get quantum corrections, and if the prediction of [29] 
is correct it would imply that these corrections drastically change the structure of 
the strong coupling expansion. 

In this paper we reconsider the computation of the Konishi state energy. We 
solve numerically the excited states TBA equations [28] for the Konishi operator up 
to 't Hooft's coupling A ~ 2046, and use the data obtained to analyze the behavior of 
Y-functions. We use the analysis to address the issue of the existence of the critical 
values of A raised in [28]. The consideration in [28] was based on an assumption that 
the analytic properties of the exact Y-functions would follow those of the large J 
asymptotic solution, and if this assumption is not realized then the TBA equations of 
[29, 28] formulated for Konishi-like states at weak coupling may be valid for all values 
of A. Extrapolating our results we find that the first critical value for the Konishi 
operator is most probably greater than 5300, A^'' > 5300, which is significantly higher 
than the estimate based on the large J asymptotic solution that gives ^ 774. 

We also find that the contribution of Yg-functions to the Konishi state energy 
grows almost linearly with the string tension g = ^/X/27r. If this behavior continues 
to hold for all values of A then this would imply that at large A the exact Bethe 
root asymptotes to a constant less than 2. This would be drastically different from 
the asymptotic behavior of the corresponding solution of the Bethe- Yang equations 
and would, in particular, imply the absence of the critical values for the Konishi 
operator. This would be a very puzzling scenario because the full spectrum of string 
theory in fiat space can be reproduced already from the Bethe- Yang equations [59], 
and if w does not asymptote to 2 it would be necessary to explain how the spectrum 
follows from the TBA equations. In particular, the spectrum degeneracy would be 
more difficult to explain because Bethe roots of different states would not behave 
uniformly at strong coupling. 

Examining the data obtained, we found convincing evidence in favor of the pre- 
diction of [29] that the first nonvanishing subleading coefficient Ci of the 1 / term 
is equal to 2, and the coefficient C2 of the terms 1/a/A vanishes. If one cuts the 
asymptotic series at order l/A^^^ and fits our data then the coefficient ci appears 
to be 2.02 ± 0.02 depending on the fitting interval used. It is clearly reasonable to 
assume that Ci is equal to 2. Setting ci = 0, one then finds that C2 = —0.02 ±0.01 
which is obviously very small. 

Then, assuming from the very beginning that c_i = Ci = 2 and cq = C2 = and 
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fitting the data with A > 77, one gets the following fitting function for the Konishi 
state energy 



/^^ oVT 2 3.26 2.53 4.03 



which differs from (1.1) by the presence of the 1/A term. Thus, the coefficient C4 of 
the term 1/A does not vanish, and we cannot confirm the prediction of [29] that up 
to an overall factor of the large A expansion of the Konishi state energy is a series 
in 1/VA. The vanishing of C2 could be related to the high degree of supersymmetry 
of the model as was pointed out in [60]. 

If on the other hand one follows [29] and assumes from the very beginning that 
C2 = C4 = then fitting the data with A > 77, one gets the following fitting function 
for the Konishi state energy 

— ifT f 1-992 2.73 7.45\ , , 

C2 = C4 = =^ Ek{\) = <f\ (^2.00005 + — + , (1.3) 

which obviously is in a very good agreement with (1.1). The coefficients in (1.3) 
mildly depend on the fitting interval used and, for example, the last coefficient can 
take values from 6 to 8. 

Let us also mention that our numerical results agree with those of [29] with the 
0.0015 precision for most values of A, and this implies the equivalence of the TBA 
equations of [11] and [28] for Konishi-like states at weak coupling.^ 

We used in our computation a new integral representation for the BES dressing 
phase [62] which significantly reduces the computational time. 

The paper is organized as follows. In the next section we present the results of 
the numerical solution of the TBA equations for the Konishi state energy. In section 
3 we discuss the properties of Y-functions and estimate the first critical value of A. In 
section 4 we use our data to find the coefficients of the large A asymptotic expansion. 
In the appendix we describe the numerical algorithm used for the computation, and 
present new formulae for various dressing phases and kernels. 



2. Konishi state energy 

We solve numerically the following equations from [28] : the simplified TBA equations 
(4.2-4.5) for Ym\wi Ym\vwi Kt -functions, and the hybrid equations (4.13) for Yq- 
functions, and determine the values of the Bethe root U2 = —Ui = w or, equivalently, 
the momentum p = p{w) carried by a string particle from the exact Bethe equation 

^Since the considerations in [28] and in v. 3 of [11] have the same starting point - the mirror 
theory string hypothesis [8] - the equivalence of the TBA equations in fact follows from the equality 
of the mirror theory dressing phases used in [29] and [28] which was recently proven in [32]. 
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(8.61). The energy of the Konishi state is then given by the following formula 



EKi\) = 2 + 2Jl + Ag^sm'^-J2 [ ^ log(l + 1q) , 



f2.11 



where p is the momentum of a Q-particle of the mirror theory, and g is the string 
tension related to 't Hooft's coupling A as A = An'^g'^. The energy is given by 
the difference of the two terms - the first one is the contribution coming from the 
dispersion relation, and the second one is the F-functions contribution. We denote 
the contributions as Edis and Ey, respectively, so that Ex = Edis — Ey. It is clear that 
Edis and —Ey play the roles of the kinetic and potential energy of the Konishi state 
particles. We will also use the notation Ey^ for the contribution of a Ig-function to 
the energy. 

It is known that at large values of 't Hooft's coupling the energy of the Konishi 
state can be expanded in an asymptotic series in powers of 1/ a/A 

E..(A) = 2^ + c„ + |L + ^ + J|_ + ^ + J|_ + .... (2.2) 

where the leading behavior follows from the string spectrum in flat space [58], 
and can be reproduced from the asymptotic Bethe ansatz [59]. The coefficient Cq is 
believed to be equal to because that is what one gets from both the free fermion 
model describing the su(l|l) sector and the asymptotic Bethe ansatz [61] but to our 
knowledge an honest string theory derivation of cq is absent. The large A perturbative 
expansion of the light-cone string sigma model (see [63] for a review) allows one to 
have any of the coefficients q nonvanishing. It was however argued in [60] that 
the coefficient C2 should vanish due to the high degree of supersymmetry of the 
model. Moreover, it follows from the fitting function (1.1) for the Konishi state 
energy obtained by solving the canonical TBA equations [29] that in fact both the 
coefficients C2 and C4 vanish, and then it is tempting to speculate that all coefficients 
with even indices vanish: C2k = 0. 

One of our aims is to demonstrate that the excited state TBA equations of [28] 
indeed predict that the leading coefficient is 2, and then to understand if one can 
fix the coefficients Ck, and in particular if the numerical data indeed predicts that 

Co = C2 = C4 = 0. 

The results of our computation of the Konishi state energy or, equivalently, the 
conformal dimension of the Konishi operator as a function of g are collected in the 
table (6.1) from the Appendix. In Figure 1 we plot the data together with the graph of 
the function -E'dis(w^asym) where iiJasym = w^asym(A) is the corresponding solution of the 
Bethe- Yang equation, and the graph of 2\/X which is the large A asymptote of both 
the exact and asymptotic energies. One sees that the exact energy graph approaches 
2 a/A faster than -E'dis(wasym)- It is also useful to plot the graph of -E'dis('U^) as a function 
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Figure 1: Black dots represent the numerical solution of the TBA equations for the Konishi state 
energy Eji{X). The brown (upper) curve represents the solution of the Bethe-Yang equation, and 
the red (lower) curve is the graph of which is the large A asymptote of Ek{^)- The range of 
the coupling constant is from g = 0.1, A ~ 0.39 to g = 7.2, A = 2046.56. 




Figure 2: Here black dots, the brown and red curves are the same as in Figure 1, and the blue 
(upper) curve is the graph of £'dis(w) where w = w{g) is the solution of the exact Bethe equations. 

of A to compare its contribution with the exact energy. The corresponding plots are 
shown in Figure 2 where we use the string tension g = \/X/27t as an independent 
variable. We see that Ed[s{w) grows almost linearly with g for g > 4: while the total 
energy grows as ^/g. This implies immediately that the Y-functions contribution 
also grows linearly, see Figure 3, and, moreover, the linear parts of E^isiw) and Ey 
cancel each other leading to the ^/g behavior of Ek- On the other hand, -Edis('Wasym) 
grows only as y/g for these values of g. It is interesting that the linear dependence 
of Ey becomes clearly visible already at very small values of g. To understand the 
reason for such a different behavior of E^j^is^w) and -Edis('U^asym) we plot in Figure 4 
the solutions of the exact Bethe equation and the Bethe-Yang equation for w and 
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Figure 4: The black dots represent our numerical solution of the exact Bethe equation, and the 
brown curve is the graph of the corresponding solution of the Bethe- Yang equation. The exact 
Bethe root w{g) reaches its minimum at g w 1, and g w 1.6 is the inflection point. 

Wasym, respectively. The numerical data of the computation of the Bethe root w are 
in the table (6.2) in Appendix. 

One sees that w is still pretty far from 2 which is the large g asymptotic value of Wasym, 
and that the exact Bethe root is noticably smaller than Wasym- The corresponding 
exact momentum p is not small at these values of g and decreases much slower than 
Pasym- This explains why the dispersion relation contribution indeed grows as g. 

It is tempting to conclude from Figures 2 and 3 that Edis{w) and Ey would grow 
linearly for all values of g, and, therefore, w in fact would asymptote to a constant 
Woo less than 2. If this happens then this would make the existence of critical values 
of the coupling constant discussed in [28] rather improbable, and as a result the TBA 



-7- 




Figure 5: The graphs of the derivative of i^dis and Ey with respect to g. 

equations for Konishi-like states [29, 28] might be valid for any value of g. Also, this 
would mean that the strong coupling limit of multi-particle states with finite number 
of particles and J (short operators in dual field theory) is very different from the near 
fiat space limit discussed in [59, 64] where the rapidities asymptote to 2, and one 
can study states with J ~ The puzzle then is that the full spectrum of string 
theory in fiat space can be reproduced already in the near fiat space limit [59], and if 
w does not asymptote to 2 it would be necessary to explain how the fiat space string 
spectrum follows from the TBA equations.^ 

This is an intriguing scenario, and it would be very interesting to understand 
analytically if it is the one. Our numerics however seems to indicate that the linear 
growth of Ey might be a feature of the intermediate coupling regime we are studying, 
and it will slow down for larger values oi g. In Figure 5 we plot the graphs of the 
derivative of E^^siw) and Ey with respect to g (obtained by using the Interpolation 
function in Mathematica). We see that the rate of change of Ey reaches its maxi- 
mum at ~ 3, remains almost constant till ~ 4 and then begins to decrease very 
slowly. We are not sure if this effect is genuine. The precision of our computation 
falls down for > 4, and the decrease in the rate of change of Ey may be just a 
numerical artifact. Since, as will be discussed in the next section, the contribution 
of an individual YQ-function to Ey slightly decreases at large g it might be also 
necessary to include the contribution of more Ig-functions than we did. Then, it is 
certainly possible that the rate of change of Ey would stabilize at even larger values 
of g, and the scenario discussed above would be realized. The assumption that the 
Y-functions contribution would be linear in g for very large g has a profound con- 
sequence on the strong coupling dependence of the Yi-function. Since Yg-functions 
are very small for |m| > 2 the integration region in (2.1) is effectively of order for 
large g, and, therefore, the Y-functions contribution would grow as g only if log Yi 
would be of order y^. However, in the next section we will see that it is not the 

^Strictly speaking, even if w asymptotes to 2 but with a rate different from the one of ifasym it 
would be a challenge to derive the flat space string spectrum from TBA. 
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ri{u,g); g= 1.6,3.0,4.4,5.8,7.2 
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Figure 6: On the left figure graphs of five Yi-functions witli different values of g are shown, and 
on the right figure the graph of Yi{0,g) as a function of g is plotted. 

case and for ~ 7, Yi increases only as g^^^. This behavior is different from both 
the asymptotic Fi-function (^-dependence and the exponential growth required by 
the scenario discussed above. This also shows clearly that the values of g we have 
reached are not large enough, and we are still analyzing the intermediate coupling 
regime. What happens at larger values of g remains to be understood. 

3. Y-functions 

In this section we discuss various properties of Y-functions. We begin with Yq- 
functions because the energy of the Konishi state depends explicitly on them. 

In Figure 6 we show plots^ of several exact Yi-functions computed at various 
values of g. One sees that Yi is getting larger in the interval [—2,2] with g increas- 
ing. In fact it increases very fast, and becomes much larger (one order of magnitude) 
than the asymptotic Y]°- function computed at the same values of g and w. In par- 
ticular, Yi{0,g) keeps increasing while the asymptotic Y]°-function computed with 
the exact Bethe root w decreases at u = for g > A. To find the (^-dependence 
we plot y/gYi{0,g) and y/gYi{l,g) in Figure 7. We see that they are almost linear 
functions, and, therefore, Yi ~ g^^'^. It is hardly possible that Yi would show the 
^f^/^-dependence for large g because this would lead to the contribution of the order 
^yglogg to the Konishi state energy which cannot be canceled by any reasonable 
contribution from E^is- If the exact Bethe root w approaches 2 for large g then 
Yi{u,g) should asymptote to a finite function. The existence of the critical values of 
g seems to require in addition that Yi{u, g) would go to for any u < 2 (but it could 
stay finite for m ~ 2 — ///(?). If w approaches Woo < 2 then, as was mentioned in the 
previous section, logFi must grow as ^/g. 

In Figure 8 the plot of the contribution i?Yi of Yi-function to the Konishi state 
energy and the plot of its derivative dgEy^ are shown. Ey^ decreases too fast to be 

^AU Y-functions are even so we plot them only for u > 0. 
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a/7 5,Fi(0,g) yfi d,Yi(l,g) 




Figure 7: The graphs of y/gdgYi{0, g) and ^/gdgYi{l, g) definitely show that Yi grows as g^^^. 




Figure 8: These are the plots of Ey^ and dgEy^- 

explained by insufficient numerical precision. The scenario with Woo < 2 would be 
realized only if dgEy^ asymptotes to a positive constant. 

Figure 9 shows plots of y2-functions. One sees that they exhibit a rather intricate 
g dependence. The maximum value of I2 increases with g and shifts to the right 
towards u = 2. 



F2(u,g); g = 1.6, 3.0, 4.4, 5.8, 7.2 i'2(0,g) 




Figure 9: On the left figure graphs of five y2-functions with difi'erent values of g are shown, and 
on the right figure the graph of 12(0,5) as a function of g is plotted. 



We see from Figure 9 that ¥2(0, g) is decreasing for g > 3. This behavior agrees 
with our expectations based on the analysis of asymptotic Y-functions. Let us recall 
that y2-function is one of the Y-functions that can be used to find the first sub- 
critical value of g because it vanishes at u = if g = g^^ [28]. Let us now assume 
that the critical value does exist. Then since Y2 is non- negative for all values of g it 
should have an expansion of the form 

Y2{Q,g)-{g--g^^f + --- . 

Since most of our points are far from gcr\ to estimate the first sub-critical value we 
first however use the linear fit. This would give the lowest estimate^ of the value of 
cjcT^ because ¥2(0, g) has a double zero at g = cjc^^ . Fitting our data in the interval 
[go, 7.2] to the function ci{g — cJct), we get the results shown in (3.1). 



go 


Fit 


50 


Fit 


50 


Fit 


6.2 
6.5 

6.8 


-0.0061(3-9.19) 
-0.0059(3-9.27) 
-0.0057(3-9.35) 


6.3 
6.6 
6.9 


-0.0060(3-9.22) 
-0.0058(3 - 9.30) 
-0.0056(3 - 9.38) 


6.4 
6.7 
7.0 


-0.0059(3 - 9.24) 
-0.0057(3 - 9.32) 
-0.0055(3 - 9.40) 



(3.1) 



As expected, the estimated value of gcr increases with go approaching 7.2, and one 
concludes from the table that gil^ > 9.4. If on the other hand one fits the data to 
C2{g - (jcr)'^, one gets 



5o 


Fit 


5o 


Fit 


5o 


Fit 


6.2 
6.5 
6.8 


0.00061(3- 11.66)2 
0.00061(3- 11.68)2 
0.00060(3- 11.70)2 


6.3 
6.6 
6.9 


0.00061(3- 11.67)2 
0.00061(3 - 11.68)2 
0.00060(3 - 11.70)2 


6.4 
6.7 
7. 


0.00061(3- 11.67)2 
0.00060(3- 11.69)2 
0.00060(3 - 11.71)2 



(3.2) 



This fitting is much more stable then the linear one, and gives cjcr ~ 11.7. It is 
certainly possible that decreasing ¥2(0, g) would slow down for larger values of g 
resulting in a larger estimate of gil\ 

The plots of Ey^ and dgEy^ are shown in Figure 10. The rate of change of 
is still increasing, and one cannot make any reliable prediction about its strong 
coupling behavior. 

Figure 10 shows similar plots of Ys-function. It is still increasing at u = but it 
is clear that it will reach its maximum soon. 



^We assume there will be no sharp changes in the behavior of 12(0,3) for larger values of 3. 




Figure 10: These are the plots of Ey^ and dgEy^. 




Figure 11: On the left figure graphs of five Ys-functions with different values of g are shown, and 
on the right figure the graph of ^3(0, 5) as a function of g is plotted. 



Even though Y3 is so small its contribution to the energy is also increasing linearly 
with g, see Figure 12, and moreover the rate of Ey-j has already reached its maximum. 




g 



1234567'=' 1234567 



Figure 12: These are the plots of Ey^ and dgEy^. 



Next we discuss Kt-functions, see Figures 13 and 14. The first subcritical value 
can be also determined from Y± because their values vanish at m = 0, and they are 
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F_(u,g); g = 1.6, 3.0, 4.4, 5.8, 7.2 y_(o,g) 




Figure 13: y_-functions 
assumed to have the following expansions^ 

Using again the linear fit, we get the following results for Y^{0,g) 



9o 


Fit 


5o 


Fit 


5o 


Fit 


6.2 
6.5 
6.8 


0.0035(3-8.93) 
0.0033(3-9.08) 
0.0030(3-9.23) 


6.3 
6.6 
6.9 


0.0034(3- 8.98) 
0.0032(3- 9.13) 
0.0030(3- 9.28) 


6.4 
6.7 
7. 


0.0033(3- 9.03) 
0.0031(3- 9.18) 
0.0029(3- 9.33) 



(3.3) 



The results in the table (3.3) are obviously compatible with those in table (3.1) 
but the estimated subcritical values gil^ in (3.3) appear to be slightly less than 
the ones from (3.1). A better estimate of gil^ is obtained by fitting the data to 



90 


Fit 


go 


Fit 


6.2 
6.4 
6.6 
6.8 


0.000045(3 - 11.14)3 ^ 0.00087(3 ^ 11-14) 
0.000042(3 - 11.25)3 + 0.00083(3 - 11-25) 
0.000040(3 - 11.36)3 + 0.00079(3 - 11-36) 
0.000038(3 - 11-46)3 + 0.00076(3 - 11-46) 


6.3 
6.5 
6.7 
6.9 


0.000043(3 - 11.19)3 + 0.00085(3 - 11-19) 
0.000041(3 - 11.30)3 + 0.00081(3 - 11.30) 
0.000039(3 - 11-41)3 ^ 0.00078(3 - 11-41) 
0.000037(3 - 11-51)3 + 0.00074(3 - 11.51) 



with the results similar to those from (3.2). 

Since l+(0, g) is still pretty far from for the values of g we are dealing with, its 
linear extrapolation to larger values of g would not give very reliable results. Indeed, 

fitting our data for Y+{0, g) in the interval [go, 7.2] to the linear function ci{g — gcr), 
we get the results shown in (3.5). 



50 


Fit 


50 


Fit 


50 


Fit 


6.2 
6.5 
6.8 


0.25(3- 9.83) 
0.24(3- 9.97) 
0.23(3-10.10) 


6.3 
6.6 
6.9 


0.25(3 - 9.88) 
0.24(3-10.01) 
0.23(3-10.14) 


6.4 
6.7 
7. 


0.24(3-9.92) 
0.23(3- 10.05) 
0.22(3 - 10.18) 



(3.5) 



Let us also mention that if one uses cs{g — (jcr Y as a fitting function one gets the 
largest of all the estimates: g^ > 16. 

Next we discuss li|^^-function, see Figure 15. This is the function that deter- 

seems possible that the expansion could be of the form Y±[0,g) ^ {g — gcr^)"^ + - - - . 
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F+(u,g); g = 1.6, 3.0, 4.4, 5.8, 7.2 




Figure 14: y+- and F_-functions are approaching at u = with g increasing. For g > gil^ they 
become positive at u = 0, and asymptote to their ground state value Y^^^*" = 1 for very large g. 



Fi|vw(u,g); g = 1.6, 3.0, 4.4, 5.8, 7.2 i'i|vw(0,g) 




Figure 15: Yi|^,„,-function is approaching — 1 at u = with g increasing. It also approaches its 
ground state value at u = oo faster for larger g. 

mines the first critical value of g because its value at m = has the following behavior 
in the vicinity of g = g^ 

Y,\,^{{),g) + l^{g-g^^f + --- . 

Using the quadratic fitting function C2{g — (jcr)'^, we get the following results 



50 


Fit 


.90 


Fit 


90 


Fit 


6.2 
6.5 
6.8 


0.0028(3-11.52)2 
0.0028(3-11.57)2 
0.0027(3 - 11.62)2 


6.3 
6.6 
6.9 


0.0028(3- 11.54)2 
0.0028(3 - 11.59)2 
0.0027(3 - 11.64)2 


6.4 
6.7 
7. 


0.0028(3- 11.55)2 
0.0027(3- 11.60)2 
0.0027(3 - 11.65)2 



(3.6) 



The results in the table (3.6) are compatible with those in table (3.2). For all values 
of go the estimated critical value gir^ appears to be less than the corresponding 
subcritical one from (3.2). They still do not differ much, and it is what one gets from 
the analysis of asymptotic Y-functions [28]. 

We conclude from the table (3.6) that gcl-^ > 11.6, and it seems reasonable 
to expect that gil^ would not exceed 12.0 (it might appear to be a too optimistic 
expectation), so the first critical value of A would be in the interval 5300 < Acr^ < 
5700. 
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Figure 17: On the left figure graphs of five Yi|u,-functions with different values of g are shown, 
and on the right figure the graph of ^11^,(0, g) as a function of g is plotted. 

The first subcritical value can be also determined from Y2\vw, see Figure 16, 
because its value at m = vanishes, and it has the following expansion 

Y2\vwi0,g) ~ g-g[l) H . 

Y2\vw{0,g) is also pretty far from and its linear extrapolation to larger values of g 
gives results slightly lower than those from table (3.2) 



.90 


Fit 


9o 


Fit 


9o 


Fit 


6.2 
6.5 
6.8 


-0.137(5-10.91) 
-0.134(5- 11.00) 
-0.131(5 - 11.08) 


6.3 
6.6 
6.9 


-0.136(5- 10.94) 
-0.133(5 - 11.02) 
-0.130(5 - 11.10) 


6.4 
6.7 
7. 


-0.135(5- 10.97) 
-0.132(5 - 11.05) 
-0.130(5-11.13) 



(3.7) 



The estimates can be made closer to the previous ones if one fits the data to higher 
order polynomials. 

Let us finally mention that y^j-functions do not show any particular gf-dependence. 
They just are increasing very fast, see Figure 17. 

It is worth stressing that the estimate of gil^ was made under an assumption that 
the first critical value exists. If it does not then Yi\^iu{0,g) would (or not) approach 
— 1 exponentially slow at (7 — t- oo. 
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4. Large A expansion from the numerical data 



As was discussed in section 2 at large values of 't Hooft's coupling the Konishi state 
energy admits an asymptotic expansion in powers of 1/ 



Ek{\) = c_r 



/~r C2 C3 



C4 C5 

+ — H — 



+ 



(4.1) 



In this section we try to understand to what extent our numerical data can be used 
to fix the coefficients Cj. We should point out however that in general an asymptotic 
series cannot be found reliably from numerical data. For example, the function 

l_eioooooo-^ 



^ v^l + A obviously asymptotes to 2\/X but any numerical computation 



l_|_elOOOOOO-yA 

performed for A < 1000000 would predict that it asymptotes to — 2v^. Thus, we have 

to assume first of all that exponentially suppressed terms become very small already 

at the values of A we are dealing with. Then, a function can approach its asymptotic 

series monotonically or in oscillations, and it does not seems possible to single out 

one from numerics. In fact, using the standard least-square fitting procedure would 

always lead to an oscillating behavior of numerical data about a fitting function. 

Next, if A is not large enough then one may need to make an assumption about the 

structure of the large A expansion, for example to decide if the series contains all 

possible terms or some of them vanish. Finally, fitting numerical data one should 

decide how many terms one should keep in an asymptotic series, and what fitting 

interval one should use. 

Since the precision of our computation is about 10""^ for (7 ~ 7 it seems reasonable 

to keep only the terms up to the 1/A^/^ order in the asymptotic expansion (4.1) . The 

fitting is done by using the data with the string tension taking values in the interval 

bo, 5^1] where gi changes from 4.0 to 7.2 (with the step 0.1), and Mathematica's Fit 

(or FindFit) functions. The first point of the fitting interval is chosen to be go = 1.4, 

A ~ 77 because it is larger than the inflection point of the Konishi state energy which 

is ^infi ~ 0.8, A ^ 25. 

We begin the fitting by making no assumption about the structure of the large 
A expansion. Below we present the table (4.2) where we use the function in (4.1) to 
fit our numerical results in (6.1). 



9i 



4. 

4.4 

4.8 

5.2 

5.6 

6. 

6.4 

6.8 

7.2 



Ai 



632. 
764. 
910. 
1070. 
1240. 
1420. 
1620. 
1830. 
2050. 



Fit 



0.0748723 
1.3609 - 
0.94197- 
0.443739 - 

-0.0777338 
-0.439602 - 
-0.131008 

0.201473 

0.210769 



19.5546 
176.565 _ 

117.655 _ 

AV4 

41.3538 _ 
A=/4 
45.6551 

A5/4 

, 110.23 , 

+ A5/4 + 

I 51.482 I 

15.5603 
AV4 
17.5595 



AV4 



_ 26.8658 
A3/4 
198.245 I 
A-V4 ^ 
137.196 

A3/4 

_ 60.5348 
A3/4 
24.225 
" A3/4 
I 85.6348 

A3/4 
I 31.0414 

A3/4 
30.0019 I 

A3/4 

31.7814 _ 
A3/4 



- 1.99888</A- 
1.94263</A- 
1.96036 v^- 

- 1.98102 v^- 
1 +2.0022a/A 

-2.01667vO^- 
f 2.00453^/1 
1. 99163 + 
1.99128 v^A- 



0.886704 
11.2876 
7.19535 

2.23404 I 
, 3.06072 
, 6.79178 
I 3.56246 
0.0366606 
0.0634158 



7.23013 
68.3193 



40.4154 



+ 



VA 
47.1638 

21.0468 , 

'^^^ 
_ 7.33942 

_ 27.6317 

\A 
__ 9.8231 _ 

/A 
I 9.86239 

10.429 

Ta 



A 

295.295 
A 

202.045 
A 

83.0606 
A 

_ 50.6135 
A 

148.678 
A 

60.4474 
A 

, 39.2596 
^ A 
42.2007 



(4.2) 
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Figure 18: On the left and right figures the graphs of c_i and Cq as functions of gi are plotted. 
The fitting is done without any constraint imposed on c^. 



To better visualize the results in Figure 18 we also plot the coefficients c_i and Cq. 
As one can see, the coefficient c_i of the leading term is oscillating about 2, and for 
gi > 6.2 it is getting quite close to 2. Its average value is 1.983. The coefficient cq is 
oscillating about but it is not really small with the average equal to 0.385. There 
is however an obvious correlation between the values of c_i and Cq - the closer c_i 
is to 2 the closer cq is to 0. 

The subleading coefficients are however not fixed at all and take very different 
values depending on the fitting interval used. The fitting therefore is not stable, and 
the fitting function strongly depends on the fitting interval. One may conclude that 
the numerical data allows one to fix only the leading coefficient in the strong coupling 
expansion. One should remember however that the series (4.1) is only asymptotic, 
and fixing its coefficients by using the data for these not very large values of A is not 
that straightforward as it is for a convergent series. 

To proceed let us fix the leading coefficient to be 2. Then, fitting the numerical 
data, one finds the results in table (4.3) and Figure 19. 
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4. 
4.4 
4.8 
5.2 
5.6 
6. 
6.4 
6.8 
7.2 



Ai 



Fit 



632. 
764. 
910. 
1070. 
1240. 
1420. 
1620. 
1830. 
2050. 



0.0487754 - 
-0.00945538 
-0.0269174- 

-0.030125 4 
-0.0216031 - 
-0.00723624 

-0.0115075 

-0.0230903 
-0.0269859 - 



16.0827 _ 
I 26.2666 
40.3386 I 

A5/4 ^ 

43.0799 I 

AV4 + 

34.8027 I 

A5/4 ^ 
I 20.0648 

1- A-V4 

I 24.7422 

AV4 
I 38.1119 
^ AV4 
42.8933 , 

' AV4 H 



23.1977 I 

A3/4 
I 7.84673 
A3/-1 
17.7895 _ 

A3/4 
19.6865 I 
A3/4 ^ 
14.2113 _ 

A3/4 

I 4.64678 

A3/1 
I 7.61427 

A3/4 
I 15.9437 
■ A-V4 
18.8607 _ 

\3/4 



2^ + 
+ 2^ 

-2^A- 
2^4 
-2^A- 
f 2^ 
+ 2^ 
f2^ 
-2v/A- 



1.13832 , 

^^"^ 
I 2.2426 

2.58159 

n 

2.64471 _ 
2.47195 

n 

2.17693 
I 2.266 
2.5106 
2.59408 



+ 



+ 



5.94492 
_ 2.36707 

Va 

_ 4.97525 

Va 

5.46707 _ 
_ 4.08328 

Va 

_ 1.69214 

Va 

2.42443 _ 

Va 

4.45848 _ 

Va 

_ 5.16215 
VA 



34.8688 
A 

22.6721 
A 

41.4559 
A 

45.0787 
A 

34.3737 
A 

15.4861 



A 

21.4161 
A 

38.2204 
A 

44.1706 



(4.3) 



The coefficient cq is now much smaller but it stops oscillating about 0. The average 
value of Co is —0.016, and Cq is decreasing with gi increasing. The subleading coeffi- 
cient Ci is not close to 2, it has an average 2.37, and is greater than 2 for almost all 
values of gi. There is still a correlation between the values of cq and ci - the closer 



-17- 



Cl(^l) 




Figure 19: On the left and right figures the graphs of cq and ci as functions of gi are plotted. 
The fitting is done with c_i ~ 2. 




Figure 20: On the left and right figures the graphs of ci and C2 as functions of gi are plotted. 
The fitting is done with c_i = 2 and cq — 0. 

Cq is to the closer ci is to 2. One sees that for the range of A we are analyzing 

the contribution of the constant term is much smaller than the one of the subleading 

term Ci/ \f\ and, therefore, it is reasonable to assume that cq = 0. 

Fixing now the leading and constant coefficients to be 2 and 0, respectively, one 
finds the results in table (4.4) and Figure 20. 
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Ai 


4. 


632. 


4.4 


764. 


4.8 


910. 


5.2 


1070 


5.6 


1240 


6. 


1420 


6.4 


1620 


6.8 


1830 


7.2 


2050 



22.2908 


1 4.1118 


A5/4 


A3/4 


18.0024 


2.18271 




A3/4 


14.458 1 


0.641502 




A3/4 


11.4777 


0.621309 


A5/4 


A3/4 


10.2483 


1.13316 


A5/4 


A3/4 


11.2064 


0.748872 


AS/4 


A3/4 


9.64851 


1.36445 


A = /4 


A3/4 


5.80801 


2.85568 


AS/4 


A3/4 


2.78511 


4.00844 


AS/4 


A3/4 



Fit 

2^fX + 
2^^X + 



2^/A+^^- 

V A 

2^- 



2.03799 

n 

2.02384 
2.01824 



2</A+^^ 

1.98967 



1.25263 _ 
0.905321 _ 
0.633195 

_ 0.413548 
0.325432 

_ 0.390144 
0.287401 
0.041104 

yx 

0.147219 



16.5214 
A 

11.8034 
A 

^ 7.96653 
A 

4.77967 



A 

3.47594 
A 

4.47424 



A 

2.86239 
A 

1.07875 
A 

4.15477 



(4.4) 



It is clear from the table and Figure 20 that fixing c_i = 2 and Cq = 
first nontrivial subleading coefficient ci to be very close to 2. Its avera; 
2.026. The next coefficient C2 still varies significantly and its average 



makes the 
ge value is 
is -0.435. 
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Figure 21: On the left and right figures the graphs of C2 and C3 as functions of gi are plotted. 
The fitting is done with c_i = 2, cq = and Ci = 2. 



However C2 is increasing and becomes very close to for = 6.9. 

Since the coefficient Ci is so close to 2, let us proceed by fixing c_i 
and ci = 2. Then, one obtains the following table and Figure 21 



2, Co 



91 


Ai 


4. 


632. 


4.4 


764. 


4.8 


910. 


5.2 


1070. 


5.6 


1240. 


6. 


1420. 


6.4 


1620. 


6.8 


1830. 


7.2 


2050. 



Fit 



6.48139 


2.7364 




la3/4 


5.92711 


2.86533 




la3/i 


5.63528 


2.93171 


la=/4 


la3/i 


5.57229 


2.94569 




la3/'' 


5.45878 


2.97011 


la=/4 


la^/* 


5.03661 


3.0599 


la'S/'i 


la^/-* 


5.00365 


3.06689 


la=/'' 


laa/'' 


5.51203 


2.96201 


la^/'' 




6.08044 


2.84643 


la=/4 





+ 2^ + 


2 


0.044888 1 


0.556775 
la 


+ 2v^4 


_ 2 


0.0331335 


1.02244 


' 


\/Ia 


la 


+ 2v^4 


_ 2 

' 


0.0271635 
\/Ia 


1.26505 

la 


+ 2v^-f 


_ 2 


0.0259249 


1.31682 




^Aa. 


^ la 


+ 2v^4 


_ 2 


0.0238023 


1.40878 




\/Ia 


^ la 


-+2^k- 


via 


0.0160542 
Via 


, 1.749 
la 


+ 2^4 


2 


0.0154523 


1.77553 




Via 


la 


+ 2v^4 


2 


0.0243275 


1.37154 




Via 


^ la 


4-2^4- 


2 


0.0340167 1 


0.922834 




VIS ^ 


la 



(4.5) 



The ffiting with the subleading coefficient equal to 2 makes the coefficient C2 to be 
very small with the average —0.026. Its contribution to the energy is much smaller 
than the contribution of the next term ca/A^^^, and it is reasonable to assume that 



C2 is equal to 0. So, let us proceed by fixing c_i = 
Then, one obtains the following table and Figure 22 



2, Co = 0, ci = 2, and C2 = 0. 
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4. 632. 

4.4 764. 

4.8 910. 

5.2 1070. 

5.6 1240. 

6. 1420. 

6.4 1620. 

6.8 1830. 

7.2 2050. 



4.19637 


3.24272 


laV4 


la^/^i 


4.15138 


3.24647 


laV4 


la3/-i 


4.11145 


3.24973 


laV4 


la^/^i 


4.05719 


3.25408 


laV4 


la^/^i 


4.0154 


3.2574 


la=/'' 


la3/-4 


4.02985 


3.25628 


laV-i 


la^/i 


4.00445 


3.25825 


laV-i 


la3/-i 


3.89399 


3.26672 


laV-i 


la3/-4 


3.75839 


3.27702 


laV-i 


la3/-i 



Fit 

+ 2^/\a- 
+ 2Vhx- 
+ 2^^ 
4-2^- 
4-2^H 
+ 2^- 
4-2^- 
4-2^- 
+ 2^- 



2 


2.43269 


via 


^ la 


2 


2.45884 


Via 


la 


2 


2.48182 


Via 


la 


2 


2.51283 


VIS 


la 


^ 4 


2.5366 


VIS ^ 


la 


2 


2.52846 


VS 


la 


2 


2.54273 


VIS 


I" la 


2 


2.60457 


VIS 


I" la 


2 


2.68017 


Via 


I" la 



(4.6) 



We see an interesting effect of this ffiting. It appears to be rather stable. As one 
can see from the table (4.5), the remaining three coefficients C3, C4 and C5 are not in 
fact too sensitive to the choice of gi anymore. 
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■3.245 
■3.250 



-3.260 
-3.265 
-3.270 
-3.275 




4.5 5.0 5^.5' 




Figure 22: On the left and right figures the graphs of C2 and C3 as functions of gi are plotted. 
The fitting is done with c_i = 2, cq = and Ci = 2. 



To continue we use the average value of C3, and fit the coefficients C4 and C5. 
Then, we use the average value of C4, and fit the coefficient C5. Finally, taking the 
average value of C5, we find the following fitting function 



^0 



1.4: EK{X) = 2<fX + 



3.26 2.53 

+ 



A 



4.03 



(4.7) 



Obviously, the fitting function is different from (1.1), in particular the coefficient C4 
is not small and gives a significant contribution to the energy. The coefficients Cj in 
(4.7) depend on the choice of Qq. In the table below we present fitting functions for 
1.2<^7o<l-6 



9o_ 
1.2 

1.3 

1.4 

1.5 

1.6 



2W+ 
2^ + 
2v^ + 
2^ + 
2^ + 



2 

2 
2 

2 
2 



Ek 


3.18 




1.89 


+ 


5.27 


A3/4 


+ 


A 


A5/4 


3.22 


+ 


2.25 


+ 


4.57 


A3/4 


A 


A5/4 


3.26 


+ 


2.53 


+ 


4.03 


A3/4 


A 


A5/4 


3.27 




2.63 


+ 


3.83 


A3/4 


+ 


A 


AS/4 


3.28 


+ 


2.70 


+ 


3.69 


A3/4 


A 


AS/4 



(4.8) 



One can see that only C3 is relatively stable C3 = — 3.23±0.05 but the other coefficients 
change more substantially. 

Does our fitting rule out the strong coupling expansion in powers of ad- 
vocated in [29]? In Figure 23 we compare the fitting function Ek{X) (4.7) with two 
fitting functions obtained by setting C2k = from the very beginning. Both Ek{X) 
and the function (1.3) fit the data equally well. This is not surprising because both 



functions have the same number of free fitting parameters. On the other hand, the 
function E,^ (A) 



'K v^v — ^ V /\ -r -2[7^ — 1^ + 1^ obtained by setting c_i = ci = 2 and 
C2k = significantly deviates from the data, up to 0.0002, which is at least twice 
more than the precision of our computation for g < 5. Thus, we are to conclude 
that the coefficient C4 does not vanish and the strong coupling expansion is in pow- 
ers of 1/ v^. Having said this, we should admit that if the actual precision of our 
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AE^(g) 



AE^(g) 



0.0002 



-0.0002 



-0.0004 



0.00005 



-0.00005 



.4- —•Ji. 
I • 



•• 6 ■ 7 



Figure 23: Black dots represent the difference between the computed values of the Konishi energy 
and its fitting function (4.7): AEx = — Ex- Green squares on the left picture represent a 



=GKV 



=GKV, 



similar difference with E^'^' where the fitting function S^"' (A) = + — |^ -I- |^ is 
obtained by using our data and setting c_i — ci = 2 and 00 = 02=04 — 0. On the right picture 
green squares represent the difference of Ex and the fitting function in (1.2). 



Ek-2 X"^ 

4 
3 
2 
1 




Ek-2 X"^ 
4 

3 

2 

1 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



8 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 



8 



Figure 24: Black dots represent the difference between the energy and its leading large A asymp- 
totic, Ex — 2v^. Red (upper) curve is the graph of Ex — 2v^. Green curve on the left picture 
is the graph of ^^^^^(A) — 2v^, and on the right picture it is the difference between the fitting 
function (1.3) and 2-^. The red curve works very well starting already with g = 0.3, A = 3.55. 



computation for g ~ 7.0 is not 0.0001 but about 0.0002, and if the contribution of 
exponentially suppressed terms is of order 0.0002 for (7 ~ 2, then this could explain 
the left plot on Figure 23 and make possible the vanishing of C4. 

Let us finally mention that the function Ek works unexpectedly well (and better 
than Ej^ or (1.3)) starting already with such a small value of the coupling con- 
stant a.s g = 0.3, A = 3.55, see Figure 24, that is less than the expected radius of 
convergency of the weak-coupling expansion, even though the fitting was done for 
the data with g > 1.4. This implicitly confirms that Ci = 2. Since -v^3.55 = 1.37 is 
close to 1, changing ci just by 0.1 would require an essential change of the coefficients 
C3, C4 and C5 to fit the energy data at large A but then the fitting at small A would 
be destroyed. 

In Figure 25 we compare our numerical results with those of [29]. 
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0.0015 
0.0010 
0.0005 




-0.0005 
-0.0010 
-0.0015 




Figure 25: The black dots represent the difference between the values of the Konishi energy 
we obtained and those from [29]. One sees that our results agree with those of [29] with 0.0015 
precision for almost all values of g. 



5. Conclusion 

In this paper we solved the TBA equations for the Konishi operator descendent from 
the s((2) sector proposed by Arutyunov, Suzuki and the author [28] up to 't Hooft's 
couphng A ~ 2046. At this value the iterations converge very slowly. It would be 
important to improve the numerical algorithm and approach closer to the first critical 
value of A, and go beyond it. One possible improvement would be to use Newton's 
method for solving the TBA equations as it was done for the Hubbard model in [65] . 
Solving the TBA equations numerically for large values of A is a challenging problem 
also because at large A the Konishi state energy is very close to the energy of the 
other two Konishi-like states with n = 1, see [28] for detail. It would be interesting 
to repeat the computation for these states, and for other Konishi-like states with 
larger string levels. 

Fitting the data for the Konishi state energy^ we found convincingly that the 
first nontrivial subleading coefficient ci is equal to 2 and the coefficient C2 vanishes 
in agreement with the previous prediction [29]. The coefficient C4 however does not 
vanish and therefore the strong coupling expansion seems to be not in powers of 
1 / -\/A (up to an overall v^) but in powers of 1 / v^- 

In principle increasing the values of A may change the fitting coefficients and 
invalidate the predictions for the coefficients. It is necessary to derive them by 
analytic means. 

Extrapolating the results obtained shows that the first critical value for the 
Konishi operator is Acr"* > 5300 and it is probably in the range 5300 < X^^ < 5700 
which is significantly higher than the estimate based on the large J asymptotic 

^We assumed the absence of log A-dependcnt terms in the asymptotic large A expansion. These 
terms are known to appear in the strong coupling expansion of the Bethe-Yang equations [66, 67]. 
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solution that gives Acr ~ 774. This means that even at such a large value of A we 
are still far from the strong coupling regime. This is an intermediate coupling regime 
and in fact it is the one where Ig-functions play the most important role. Moreover, 
if the contribution of Ig-functions to the Konishi state energy will continue growing 
as then this would imply that at large A the exact Bethe root asymptotes to a 
constant less than 2, and the critical values for the Konishi operator are absent. 

It is clear that the numerical computation we performed raised more questions 
than gave answers. Some of the questions can be answered only analytically, and we 
hope to address them in future. 
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6. Appendix 

Numerical data 

Here we collect our numerical data. 

In the table (6.1) we present the results of the computation of the energy of the 
Konishi state or, equivalently, the conformal dimension of the Konishi operator as a 
function of g 
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In the table (6.2) we present the results of the computation of the Bethe root w 
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Numerical algorithm 

We compute the Konishi state energy in several steps. At the first step one solves the 
TBA equations for a fixed Bethe root w by iterations. Equations for Kt -functions are 
solved first, then equations for Y^, and and finally equations for Yq. For g = 0.7 
at the first iteration one uses the vacuum solution for Yy^- and Yuto-functions [31], 
and asymptotic YQ-functions [38].^*^ For larger g the solution found at the previous 
value of g, and the linear extrapolation of the Bethe root are used. 

The number of iterations is bounded not to exceed Alitor = 10, and the iterations 
also stop if the absolute value of the difference between the energies of two successive 
iterations, dEn = \En — -E'n-il, becomes less than de/10 where de = 0.01 is the 
precision the Konishi state energy is computed with at this step. 

The number of Yg-functions computed is determined by the contribution of the 
asymptotic Yg-functions to the energy. It is given by J dp log(l + Yq), and we solve 
the TBA equations only for those Yg-functions whose contributions exceed 10~^. We 
denote the total number of Yg-functions used in the computation by Qmax? and it 
depends on the value of g, to be precise we used, Qmax = 3 for g = 0.7, Qmax = 4 for 
0.8<g< 1.1, Q^ax = 5 for 1.2 < g < 1.8, Q^^x = 6 for 1.9 < g < 2.8, Q^r.. = 7 for 
2.9 < ^ < 4.2, and Q^ax = 8 for 4.3 < ^ < 7.2. 

The number of Y^-functions is equal to = 15, and number of Y^u,-functions 
is equal to + Qmax — 2. This seems to be more than we need due to the locality of 

^"For g < 0.7 wc just use the asymptotic solution for Yg-functions. It is sufficient for the precision 
we are after. For smaU enough values of g one can also use the asymptotic large J solutions for 
Y-functions to start iterations. They are expressed in terms of transfer matrices corresponding 
to various representations of the symmetry algebra of the model under consideration [68, 69] (see 
also [30]). In the AdSs x case the symmetry algebra of the light-cone string theory is the 
centrally extended su(2|2) superalgebra [70, 71], and explicit formulae for the transfer matrices 
were conjectured in [72] and derived in [73]. 
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the simplified TBA equations for auxiliary functions, and because the hybrid TBA 
equations for Ig-function involve only Yi|„^- and Ig^ii^u^-functions. 

Y-functions are computed at discrete sets of points. In particular, YQ('u)-functions 
are computed at points Uk = kdu , k = 0,1, ... with the step du = 10/1001 in the in- 
terval [0, where u^^^{Q) is the lowest value of u such that Yqiu'^axiQ)) < 
10^^. Then, one uses Mathematica function Interpolation, and also the large u behav- 
ior of asymptotic Yg-functions to define Yq as continuous functions on the real w-line. 
Yu;- (and Yy^-functions) are computed until |1 — Yq|^(mJ^^^(Q))/((5(Q + 2))| < 10"^ 
with a variable step equal to du for < m < 4, to 2du for 4 < -u < 12 and so on. 
One uses again Interpolation and the large u asymptotic of Y^,- and F^^-functions 
to define them for all u. Finally, Kt ('u)-functions are computed at points = kduy 
with the step duy = 1/10. 

At the second step one uses the Y-functions to solve the exact Bethe equation, 
and find an adjusted value of w and energy. Then one repeats the procedure until 
the difference between the energies becomes less than de = 0.01. 

At the third step one uses the solution obtained to refine it by decreasing the 
steps du and dUy, and the energy precision de. We first solved the TBA equations 
with du ^ 1/50, duy = 1/50 and de = 0.0001. Then we decreased the steps to 
du ~ 1/100, duy = 1/100. Finally we used the steps du ~ 1/130, duy = 1/130 for 
3.9 < g < 4.9, and du ~ 1/150, duy = 1/150 for g > 5.0. The precision is about 
0.0005 for g < 5.0. For larger values of g the precision decreases and is about 0.001. 

In version 2 of the paper to increase the precision we switched from the u variable 
to the mirror momentum pQ rescaled so that the intervals [0, 2] are mapped to each 
other. We then used (^-dependent steps dp = 2/NQ^/2/g with A^q = 200 and Nq = 
250. Finally, for (7 > 3 we increased the number of auxiliary functions A^w = 20. 
Observing the changes in the energy, we believe that the final precision is about 
0.0001 for all values of g. 

To solve the hybrid TBA equations^^ we compute the dressing kernels Kq-^^{u, v) 
at the points {kdu, ndu) , k,n = 0,1, ... with the step du = 50/1001, and Kqq,{u, v), 
Q' >2 with the step du = 100/1001 in the rectangle [0,ul^^^{Q)] x [0, and 
KqiS'^^'^) at the points {kdu,Av) , k = 0,1, where Av contains the Bethe root 
w, and the step in the u-direction is dv = 1/2000. 

Dressing phases and kernels 

The improved dressing phases in various kinematic regions of the string and mirror 
models involve a common function which can be written as a sum of ^-functions, see 
[74] for definitions. We denote this function as Q{xf ,Xi,ui,Q,X2,X2,U2,M,g), 

^^To avoid computing many dressing kernels we first tried to use the simplified TBA equations 
for Yg-functions. It appeared however that they were numerically unstable (at least with our 
algorithm) and led to a large systematic error increasing with g. 
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where xf satisfy the usual constraints 

1 i 1 i 

xf + ^ = ui±-Q, xf + ^ = U2±-M. (6.3) 
x^ g X2 g 

If the parameters xf belong to the string region, that is |a;^| > 1, then O is just 

equal to the BES dressing phase. 

One can show by using the DHM representation [75] for the dressing phase that 
for any choice of xf satisfying (6.3), has an integral representation involving elliptic 
H-function. Since the derivation of the representation is rather involved we present 
below the expression leaving a proof of the formula to curious postdocs 

i {{x^Y ~ 1) (2^2" ~ 2;J)(a;^a;J(5 + Q + ig{x2X^Ui + lii - 2{2x^x^{l - t) + 2(1 - t) + x^ + x^))) 
{-iQ - 2g{l -t)+ gui) {AP - (Q + ig(-4(l - t) + ui - uz))^) x^x^x+ 

.{-iQ-2g{l-t)+gu,Y^^l~^) 
i {{x'^Y — 1) — X2){x2X2Q + Q — ig{x2X2Ui +ui - 2(2x2" 2:^(1 ^ + 2(1 -t)+X2+ X2))) 



xn 



{iQ - 2g{l -t)+ gui) {AP - (Q - ig(-4(l - t) + ui - U2W) X2xlxJ 

,{iQ-2g{l-t)+gui)^\2-t)^J 
{{x^)'^ - 1) {x2 - x^){iQ{x2xj + 1) + 5(4(1 -t) + ui- 2x2 + (4(1 - t)x2 + U1X2 - 2)x^j) 
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(Q2 _ (M - ig{4{l -t) + ui - U2)y) {-iM - 2g{l -t)+ gu2)x^X2X+ 
4g2t2 e X 
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This integral representation significantly speeds up computing the dressing phases. 
For small values of g it might be useful to do the change t ^ {1 — t)/g but for very 
large g the representation seems to be the best one for numerics. 

The function 9 is then used to define the corresponding expressions Qgs, ©sm, 
Qms and Qmm in different regions of the 2;-torus 

Qss{u,Q,v,M,g) = e (xs(u+ —),Xs{u~ —),u,Q,Xs(v + —),Xs(v- —),v,M,g] (6.5) 
V 9 9 9 9 J 

Qsm{u,Q,v,M,g) = e [xsiu+ —),Xs{u'- —),u,Q,x(v + —),x{v- —),v,M,g] 

V 9 9 9 9 J 

Q,nsiu,Q,v,M,g) = e [x(u+ —),x(u~ —),u,Q,Xs(v+ —),Xs(v- —),v,M,g] 

V 9 9 9 9 J 

Qm7n{u,Q,v,M,g) = e (x(u+ —),x(u- —),u,Q,x(v+ —),x(v'- —),v,M,g] 
\ 9 9 9 9 J 

Here x{u) and Xs{u) are the mirror and string x-functions: x{u) = |(m — i\/4 — m^), 
and Xs{u) = ^u{l + \/l — A/u^). In particular, Qgg is just the BES dressing phase. 

The improved dressing phases also involve contributions of \l/-functions [74]. The 
total contribution of ^'(x2, a;i)-functions to the improved dressing phase in the string- 
mirror region is given by 

^2i{x+,x^,U2,g,M)= desm{e) (6.6) 



igM ^ C —i{x-^ — x^) siii(6') — (xj^ + x'^) cos[0) + x-^ x^ + \ 



Tz(^A[a^-gcos{e)f + AP^ V K^i - .Tf)sin(6l) - (xj + x+) cos{e) + x ^ x+ + I 
2M ( -i{x^ -x^)s\n{e) + {x^ +x^)cos{e)+x^x^ + 1 



A{S^+gcos{e)f + AP V - x+)s\n(e) + {x^ + x+) cos{e) + x^ x+ + 1 

y + ^-,.cos(.) + l)+V.(°) (t-'?+^ 
—i{xi — x'l) su\{9) — {x^ + Xi) cos{6) + x^x^ + 1 



X lo2 



X log 



i{x^ — x^) sm{9) — {xi + x^) cos{9) + x^ xf + 1 

f -^-.9cos(.) + l)+0(o) (f + ^+. 

-i{xi — x^) sin(0) + {x^ + x^) cos(0) + x^ x^ + 1 
i{xi — xf) sm{9) + {xi + Xi) cos{9) + x^xf + 1 



Here the definitions are the same as above but Q = 1, and xj^ are in the string 
region, that is xf = Xs{ui ± ^), while xf are in the mirror region, that is xf = 

x{u2 ± g^)- It is convenient to introduce 

^'^{ui,Q,U2,M,g)^^2i [xs{ui + ^-j),Xs{ui-^-^),U2,g,M^ (6.7) 
Then, the contribution of the last term in the formula for the improved string-mirror 
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dressing phase is given by 



jlogi:Zt{ui,Q,U2,M,g)= (6.8) 

7^ log ^ ^ ^ — 

(-.(«.-f)-^) («.(«.+f)-h+f)) 

The total improved dressing phase | log Ei^m in the string-mirror region is thus 
given by the sum of the three terms 

-logSi.M(«i,l,«2,M,s() = e'''"(ui,l,ii2,M,c,) + *^\"(ui,l,U2,M,5f) 
t 

+ ^logSSt(wi,l,n2,M,ff) (6.9) 



The improved dressing phase in the mirror-string region is related to the one in 
the string-mirror region by the unitarity relation 



-logEMu{u2,ui) = --logEi^M(tii,ti2) ■ (6.10) 



Finally, the improved dressing phase in the mirror-mirror region is given by the 
sum of four terms. The first one is Qrnm- The second term is the total contribution 
of ^(a;2, xi)-functions given by 

2 r (ig - ^i(gu2 - iM) + Ij r (ig - ^i(iM + gu^) + Ij 
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The third term is the total contribution of "^{xi, X2)-functions given by 

(u. x+ X- oM)^ ' loe r i-W + ¥(9^1 - ^M) + 1) r {-^g + ^^{^M + gu,) + l) 



(6.12) 
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We also introduce 

^Tr{ui,Q,U2,M,g) = ( U2 

Finally, the last term is given by 
logI]i^-(ui,Q,U2,M,.g) = 



iM 



,X \U2 



iM 
9 



.9, Q 



(6.13) 



X Ui 



(6.14) 



log 



^QT {M - \^{^{Q M) + gu, gu2)) (l 



^ ^MT (Q + \^{^{Q - M) + gu, - .g.^)) (l - .(^^^m.^^.,,^) ) ^ (' 

Thus, the improved dressing phase Sqm in the mirror-mirror region is given by 

■ log Sqm {u,,Q,U2,M,g) = e'^'^iui,Q,U2, M, 5) + (ui , Q, ^2 , M, 5) + ("1 , Q, "2 , M, <?) 



1 



log5]S,'?(«i,Q,W2,M,5) 



(6.15) 



The formulae above are valid for any real mi, M2, and should be used with caution 
for complex values of the parameters. 

Finally, the dressing kernels are defined as 



KnAAu.v) = J—-^\0KT,om(u,v) , Kn^ (u, v) = J— logT^qi^U, v) . 



Mathematica form of the dressing phases 



To make sure there are no misprints in the formulae above we also present below the 

same expressions in Mathematica input form. 
Eq. (6.4) takes the form 
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\ [CapitalTheta] [xpl_, xml_, ul_, Q_, xp2_, xm2_, u2_, M_, g_] := 

(1/Pi"2) *I*g'2*NIiitegrate [ ( ( ( (-l+xml"2) * (xm2-xp2) ) / (xml*xm2*xp2) ) * 

( ( (-1) *q* ( l+xm2*xp2) +g* (4-4*t+ul-2*xm2+ (-2+4* ( 1-t ) *xm2+ul*xm2) *xp2) ) / 

((I*Q-g*(2*(l-t)+ul))*(M~2-(Q+I*g*(4*(l-t)+ul-u2))'2)))* 

EllipticPi[((2*g*t)/(I*Q-g*(2*(l-t)+ul)))-2, (t/(2-t))-2]+ 

(- ( ( (-l+xpl"2) * (xm2-xp2) ) / (xm2*xpl*xp2) ) ) * ( (I*Q* (l+xm2*xp2) + 

g*(4*(l-t)+ul-2*xm2+(-2+4*(l-t)*xm2+ul*xm2)*xp2))/(((-I)*Q-g*(2*(l-t)+ul))* 

(M"2-(Q-I*g*(4*(l-t)+ul-u2))~2)))* 

EllipticPi[((2*g*t)/((-I)*q-g*(2*(l-t)+ul)))-2, (t/(2-t))-2] 
+(-((!* (-l+xpl~2) * (xm2-xp2) ) / (xm2*xpl*xp2) ) ) * 

( (Q+C)*xin2*xp2-I*g* (ul+ul*xin2*xp2-2* (2* (1-t) +xm2+xp2+2* (1-t) *xm2*xp2) ) ) / 
( (I*Q-2*g* (1-t) +g*ul) * (M~2- (Q-I*g* (-4* (1-t) +ul-u2) ) ~2) ) ) * 
EllipticPi[((2*g*t)/(I*Q-2*g*(l-t)+g*ul))"2, (t/(2-t))~2]+ 
(- ( (I* (-l+xml~2) * (xm2-xp2) ) / (xml*xm2*xp2) ) ) * 

( (q+Q*xm2*xp2+I*g* (ul+ul*xm2*xp2-2* (2* (1-t) +xm2+xp2+2* (1-t) *xm2*xp2) ) ) / 

( ( (-1) *q-2*g* ( 1-t) +g*ul) * (M~2- (q+I*g* (-4* (1-t) +ul-u2) ) "2) ) ) * 

EllipticPi[((2*g*t)/((-I)*q-2*g*(l-t)+g*ul))"2, (t/(2-t))~2]+ 

(- ( ( (xml-xpl) * (-l+xp2"2) ) / (xnil*xpl*xp2) ) ) * ( (I*M* (l+xml*xpl) + 

g* (4* (1-t) +u2-2*xml+ (-2+4* (1-t) *xml+u2*xml) *xpl) ) / ( (I*M+2*g* (l-t)+g*u2) * 

(q-2-(M+I*g*(-4*(l-t)+ul-u2))-2)))* 

EllipticPi [ ( (2*g*t) / (I*M+2*g* (1-t) +g*u2) ) "2 , (t/ (2-t) ) ~2] 
+ ( ( (-l+xm2"2) * (xml-xpl) ) / (xml*xm2*xpl) ) * 

(((-I)*M*(l+xml*xpl)+g*(4*(l-t)+u2-2*xnil+(-2+4*(l-t)*xml+u2*xml)*xpl))/ 
( (q~2- (M-I*g* (-4* (1-t) +ul-u2) ) "2) * ( (-1) *M+2*g* (1-t) +g*u2) ) ) * 
EllipticPi [ ( (2*g*t) / ( (-1) *M+2*g* (1-t) +g*u2) ) "2 , (t/ (2-t) ) ~2] + 
( ( (xml-xpl) * (-l+xp2~2) ) / (xml*xpl*xp2) ) * ( (I*M* (l+xml*xpl) + 

g* (u2+u2*xml*xpl-2* (2* (1-t) +xml+xpl+2* (1-t) *xml*xpl) ) ) / ( (I*M-2*g* (1-t) +g*u2) * 
(q"2-(M+I*g*(4*(l-t)+ul-u2))~2)))* 

EllipticPi [ ( (2*g*t) / (I*M-2*g* (1-t) +g*u2) ) '2 , (t/ (2-t) ) ~2] 
+ (- ( ( (-l+xm2'2) * (xml-xpl) ) / (xml*xm2*xpl) ) ) * 

( ( (-1) *M* (l+xml*xpl) +g* (u2+u2*xml*xpl-2* (2* (1-t) +xml+xpl+2* (1-t) *xml*xpl) ) ) / 
( ( (-1) *M-2*g* ( 1-t) +g*u2) * (q~2- (M-I*g* (4* (l-t)+ul-u2) ) ~2) ) ) * 
EllipticPi [ ( (2*g*t) / ( (-1) *M-2*g* (1-t) +g*u2) ) "2 , (t/ (2-t) ) ~2] ) * 
((LogGamma[l-2*I*g*(l-t)]-LogGamma[l+2*I*g*(l-t)] )/(2-t)) ,-[t,0,l}] ; 

Here xpl = xf, and so on, xf + l/xf = Ui + -Q, and Xg" + = M2 + -M. For 
effective numerics one should exclude the point t = 1 — M2I (if ^^2! < 4) from 
the integration region. This doubles the formula above, and by this reason we leave 
the necessary modification of the formula to the interested reader. For some points 
close to ±2 the formula works only if one specifies WorkingPrecision. For small values 
of g and for the precision we are after it is sufficient to use WorkingPrecision— j-lO. 

©ss, ©sm, ©ms and Qmm in (6.5) are written as 

\ [CapitalTheta] ss [u_, q_, v_, M_, g_] := 

\ [CapitalTheta] [xs [u+(I/g)q] ,xs [u-(I/g)q] ,u,q,xs [v+(I/g)M] ,xs [v-(I/g)M] ,v,M,g] ; 
\ [CapitalTheta] sm [u_ , q_, v_, M_, g_] := 

\ [CapitalTheta] [xs [u+(I/g)q] ,xs [u-(I/g)q] ,u,q,x[v+(I/g)M] ,x[v-(I/g)M] ,v,M,g] ; 
\ [CapitalTheta] ms [u_, q_, v_, M_, g_] := 

\ [CapitalTheta] [x[u+(I/g)q] ,x[u-(I/g)q] ,u,q,xs [v+(I/g)M] ,xs [v-(I/g)M] ,v,M,g] ; 
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\ [CapitalTheta]mm[u_, Q_, v_, M_, g_] := 

\ [CapitalTheta] [x[u+(I/g)Q] ,x[u-(I/g)Q] ,u, Q ,x [v+(I/g)M] ,x[v-(I/g)M] ,v,M,g] ; 



Here x[u] and xs[u] are the mirror and string x-functions. 

The total contribution of \l/(a;2, a;i)-functions to the improved dressing phase in 
the string-mirror region in (6.6) is given by 

\ [CapitalPsi] 21 [xpl_,xinl_,u2_,g_,M_] := g/(2 \ [Pi] I) (Nlntegrate [ 
Sin[\[Theta]] (Log[(l+xinl xpl-(xinl+xpl) Cos [\ [Theta] ] 
-I (xml-xpl) Sin [\ [Theta] ])/(l+xinl xpl-(xml+xpl) 

Cos [\ [Theta] ] +1 (xml-xpl)Siii[\ [Theta] ] )] ( (-2 M)/(M-2+4((g u2)/2-g Cos [\ [Theta] ] ) "2) ) ) 
+ Sin [\ [Theta]] (Log [(1+xml xpl+(xml+xpl) Cos [\ [Theta] ] 
-I (xml-xpl) Sin [\ [Theta] ])/(l+xinl xpl+(xml+xpl) 

Cos [\ [Theta] ]+I (xml-xpl )Sin[\ [Theta] ])]( (-2 M)/(M-2+4((g u2)/2+g Cos [\ [Theta] ]) "2) ) ) 

+ Sin [\ [Theta]] Log [(1+xml xpl- (xml+xpl) Cos [\ [Theta] ] 

-I (xml-xpl) Sin [\ [Theta]]) /(1+xml xpl-(xml+xpl) 

Cos [\ [Theta] ]+I (xml-xpl) Sin [\ [Theta] ] )] 

(PolyGamma[l+M/2+(I g u2)/2-I g Cos [\ [Theta] ] ] 

+PolyGamma[l+M/2-(I g u2)/2+I g Cos [\ [Theta] ]] ) 

+Sin[\ [Theta]] Log [(1+xml xpl+(xml+xpl) Cos [\ [Theta] ] 

-I (xml-xpl) Sin [\ [Theta]]) /(1+xml xpl+(xml+xpl) 

Cos [\ [Theta] ]+I (xml-xpl) Sin [\ [Theta] ] )] 

(PolyGamma[l+M/2+ (I g u2)/2+I g Cos [\ [Theta] ] ] 

+PolyGamma[l+M/2-(I g u2)/2-I g Cos [\ [Theta] ]]), {\ [Theta] ,0 , Pi/2}]); 

Here the definitions are the same as above but Q = 1, and xf are in the string 
region, that is xf = Xs{ui ± while xf are in the mirror region, that is xf = 

x{u2 ± g^)- Then, (6.7) is written as 

\[CapitalSigma]\ [CapitalPsi] 21 [ul_, Q_, u2_, M_, g_] := 

\ [CapitalPsi] 21 [xs[ul + I*(Q/g)], xs [ul - I*(Q/g)], u2, g, M] ; 

The last term in the formula (6.8) for the improved string-mirror dressing phase 
is given by 

\[CapitalSigma]lastsm[ul_, Q_, u2_, M_, g_] = (1/ (2*1) ) *Log [(xs [ul - I*(Q/g)] 

- x[u2 + I*(M/g)])*(xs[ul - I*(q/g)] - l/x[u2 - I*(M/g)])*((xs[ul + I*(Q/g)] 

- l/x[u2 - I*(M/g)])/((xs[ul + I*(Q/g)] - x[u2 + I*(M/g)] )*(xs [ul - I*(Q/g)] 

- l/x[u2 + I*(M/g)])-2))] ; 

The total improved dressing phase (6.9) in the string-mirror region is thus given 
by the sum of the three terms 

\ [CapitalSigma] Ism[ul_, 1, u2_ , M_, g_] := \ [CapitalTheta] sm [ul , 1, u2, M, g] 

+ \ [CapitalSigma] \ [CapitalPsi] 21 [ul,l,u2,M,g] + \ [CapitalSigma] lastsm [ul , 1 ,u2 ,M,g] ; 

The total contribution (6.11) of \E'(a;2, a;i)-functions in the mirror-mirror region 

is 
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\ [CapitalPsi] 21m[xpl_, xml_, u2_, g_, M_] := 
(g/(2*Pi))*NIntegrate[Sin[\[Theta]]*(2*\[Theta] - 2*Pi 

- I*Log[(l - l/(E~(I*\[Theta])*xml))/(l - E~(I*\[Theta] )/xml)] - 

I*Log [ (1-xpl/E- (I*\ [Theta] ) ) / (1-E~ (I*\ [Theta] ) *xpl) ] ) * ( (-2*M) / (M-2+4* ( (g*u2) /2 
- g*Cos[\ [Theta] ])~2) + PolyGammaCl + M/2 + (I*g*u2)/2 - I*g*Cos [\ [Theta] ] ] 
+ PolyGammaCl + M/2 - (I*g*u2)/2 + I*g*Cos [\ [Theta] ]] ) + 

Sin [\ [Theta] ] * (-2*\ [Theta] -I*Log [(1+E~ (I*\ [Theta] ) /xml) / ( 1+1/ (E" (I*\ [Theta] ) *xml) ) ] - 
I*Log [ d+E" [Theta] ) *xpl) / (l+xpl/E~ (I*\ [Theta] ))])*( (-2*M) / (M'2+4* ( (g*u2) /2 
+g*Cos[\ [Theta] ])~2) + PolyGamma[l + M/2 + (I*g*u2)/2 + I*g*Cos [\ [Theta] ] ] + 
PolyGainma[l + M/2 - (I*g*u2)/2 - I*g*Cos [\ [Theta] ]]) , {\ [Theta] , 0, Pi/2}] + 
(1/2) * (LogGamma [1+(I* (g*u2+I*M) ) /2-I*g] -LogGamma [1- (I* (g*u2+I*M) ) /2+I*g] + 
LogGamina[l+(I*(g*u2-I*M))/2-I*g] -LogGainma[l-(I*(g*u2-I*M))/2+I*g] ) ; 

The total contribution (6.12) of \l/(a;i, a;2)-functions is given by 

\ [CapitalPsi] 12m [ul_, xp2_, xm2_, g_, M_] := 

(g/(2*Pi))*NIntegrate [Sin [\ [Theta] ]*(-2*\ [Theta] + 2*Pi 

+ I*Log[(l - l/(E~(I*\[Theta])*xm2))/(l - E~ (I*\ [Theta] )/xm2)] + 

I*Log [ (l-xp2/E~ (I*\ [Theta] ) ) / (1-E~ (I*\ [Theta] ) *xp2) ] ) * ( (-2*M) / (M"2+4* ( (g*ul) /2 

- g*Cos[\ [Theta] ])~2) + PolyGaimna[l + M/2 + (I*g*ul)/2 - I*g*Cos [\ [Theta] ] ] 

+ PolyGaiimia[l + M/2 - (I*g*ul)/2 + I*g*Cos [\ [Theta] ]] ) + Sin [\ [Theta] ]* (2*\ [Theta] 
+ I*Log[(l + E~(I*\[Theta] )/xm2)/(l + 1/ (E~(I*\ [Theta] )*xm2))] 

+ I*Log [ (1+E~ (I*\ [Theta] ) *xp2) / (l+xp2/E" (I*\ [Theta] ))])*( (-2*M) / (M~2+4* ( (g*ul) /2 
+ g*Cos[\ [Theta] ])~2) +PolyGamma[l + M/2 + (I*g*ul)/2 + I*g*Cos [\ [Theta] ] ] 
+ PolyGainma[l + M/2 - (I*g*ul)/2 - I*g*Cos [\ [Theta] ]]) , {\ [Theta] , 0, Pi/2}] - 
(1/2) * (LogGamma [1+ (I* (g*ul+I*M) ) /2-I*g] -LogGamma [1- (I* (g*ul+I*M) ) /2+I*g] + 
LogGamina[l + (I*(g*ul - I*M))/2 - I*g] -LogGamma[l - (I*(g*ul - I*M))/2 + I*g]); 

^^'^ and ^^1°^ functions in (6.13) are 

\ [CapitalSigma]\ [CapitalPsi] 12m [ul_, Q_, u2_, M_, g_] 

:= \ [CapitalPsi] 12m [ul, x[u2 + I*(M/g)], x [u2 - I*(M/g)], g, Q] ; 

\[CapitalSigma]\ [CapitalPsi] 21m [ul_, Q_, u2_, M_, g_] 

:= \ [CapitalPsi] 21m [x[ul + I*(q/g)], x[ul - I*(Q/g)], u2, g, M] ; 

Finally, the last term (6.14) is given by 

\ [CapitalSigma] lastmm [ul_ , Q_, u2_, M_, g_] := 

(l/I)Log[(I'q/I"M) (Gamma [M- (1/2) (g*ul-g*u2+I (Q-M) )] /Gamma [Q+ (1/2) (g*ul-g*u2+I(Q-M))] )* 
((l-l/(x[ul+I*(Q/g)]*x[u2-I*(M/g)]))/(l-l/(x[ul-I*(Q/g)]*x[u2+I*(M/g)])))* 
( (x [ul+I* (Q/g) ] /x [ul-I* (Q/g) ] ) * (x [u2-I* (M/g)] /x [u2+I* (M/g) ] ) ) " (1/2) ] ; 

Thus, the improved dressing phase (6.15) in the mirror-mirror region is 

\ [CapitalSigma] Imm[ul_, Q_, u2_, M_ , g_] := 
\ [Capit alThet a] mm [ul , Q , u2 , M , g] +\ [CapitalSigma] \ [CapitalPsi] 12m [ul , Q , u2 , M , g] 
+\ [CapitalSigma] \ [CapitalPsi] 21m [ul,Q,u2,M,g]+\ [CapitalSigma] lastmm [ul,Q,u2,M,g] ; 
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